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We investigate the two-word Naming Game on two-dimensional random geometric graphs. Studying this 
model advances our understanding of the spatial distribution and propagation of opinions in social 
dynamics. A main feature of this model is the spontaneous emergence of spatial structures called opinion 
domains which are geographic regions with clear boundaries within which all individuals share the same 
opinion. We provide the mean-field equation for the underlying dynamics and discuss several properties of 
the equation such as the stationary solutions and two-time-scale separation. For the evolution of the opinion 
domains we find that the opinion domain boundary propagates at a speed proportional to its curvature. 
Finally we investigate the impact of committed agents on opinion domains and find the scaling of consensus 
time. 

Relevant features of social and opinion dynamics 1 3 can be investigated by prototypical agent-based models 
such as the voter model 4,5 , the Naming Game 6,7 , or the majority model 8,9 . These models typically include a 
large number of individuals, each of which is assigned a state defined by the social opinions that it accepts 
and updates its state by interacting with its neighbors. Opinion dynamics driven by local communication on 
geographically embedded networks is of great interest to understanding the spatial distribution and propagation 
of opinions. In this paper we investigate the Naming Game (NG) on random geometric graphs as a minimum 
model of this type. We focus on the NG but will also compare it with other models of opinion dynamics. 

The NG 6,7 was originally introduced in the context of linguistics and spontaneous emergence of shared 
vocabulary among artificial agents 10,11 to demonstrate how autonomous agents can achieve global agreement 
through pair-wise communications without a central coordinator. Here, we employ a special version of the NG, 
called the two-word 1215 Listener-Only Naming Game (LO-NG) 15,16 . In this version of the NG, each agent can 
either adopt one of the two different opinions A, B, or take the neutral stand represented by their union, AB. In 
each communication, a pair of neighboring agents are randomly chosen, the first one as the speaker and the 
second one as the listener. The speaker holding A or B opinion will transmit its own opinion and the neutral 
speaker will transmit either A or B opinion with equal probability. The listener holding A or B opinion will 
become neutral when it hears an opinion different from its own and the neutral listener will adopt whatever it 
hears. Detailed instances are shown in Supplementary Table. 

Consensus formation in the original NG (and its variations) on various regular and complex networks have 
been studied 6,7,1317 22 . In particular, the spatial and temporal scaling properties have been analyzed by direct 
simulations and scaling arguments in spatially-embedded regular and random (RGG) graphs 17,20,21 . These results 
indicated 17,20,21 that the consensus formation in these systems is analogous to coarsening 23 . In this paper, we 
further elucidate on the emerging coarsening dynamics in the two-word LO-NG on RGG by developing mean- 
field (or coarse-grained) equations for the evolution of opinions. Our method of relating microscopic dynamics to 
macroscopic behavior shares similar features with those leading to effective Fokker-Planck and Langevin equa- 
tions in a large class of opinion dynamic models (including generalized voter models with intermediate states) 24,25 . 

A random geometric graph (RGG), also referred to as a spatial Poisson or Boolean graph, models spatial effects 
explicitly and therefore is of both technological and intellectual importance 26,27 . In this model, each node is 
randomly assigned geographic coordinates and two nodes are connected if the distance between them is within 
the interaction radius r. Another type of network with geographic information is the regular lattice. Fundamental 
models for opinion dynamics on regular lattice has been intensively studied 3,4,28 . In many aspects, opinion 
dynamics behaves similarly on RGGs and regular lattices with the same dimensionality, but in our study, we 

also observe several differences. For example, the length scale of spatial coarsening for large f is lit) = P /in t on 
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RGG while it is P on regular lattice. More generally, concerning the 
spatial propagation of opinions in social systems or agreement 
dynamics in networks of artificial agents, random geometric graphs 
are more realistic for a number of reasons: (/) RGG is isotropic (on 
average) while regular lattice is not; (it) the average degree (k) for an 
RGG can be set to an arbitrary positive number, instead of a small 
fixed number for the lattice; (Hi) RGGs closely capture the topology 
of random networks of short-range-connected spatially-embedded 
artificial agents, such as sensor networks. 

An important feature of the NG which makes it different from 
other models of opinion dynamics, e.g., voter model, is the spontan- 
eous emergence of clusters sharing the same opinion. Generally these 
opinion clusters are communities closely connected within the net- 
work. This feature of the NG can be used to detect communities of 
social networks 19 . For the NG on random geometric graph concerned 
in this paper, the clusters form a spatial structure to which we refer to 
as opinion domains and which are geographic regions in which all 
nodes share the same opinion. A number of relevant properties of the 
NG on random geometric graph have been studied by direct indi- 
vidual-based simulations and discussed in 20,21 , such as the scaling 
behavior of the consensus time and the distribution of opinion 
domain size. In contrast to these previous works, here we develop 
a coarse-grained approach and focus on the spatial structure of the 
two-word LO-NG, such as the correlation length, shape, and pro- 
pagation of the opinion domains. 

In this paper, we provide the mean-field (or coarse-grained) equa- 
tion for the NG dynamics on RGG. By analyzing the mean-field 
equation, we list all possible stationary solutions and find that the 
NG may get stuck in stripe-like metastable states rather than achieve 
total consensus. We find significant two-time-scale separation of the 
dynamics and retrieve the slow process governed by reaction-dif- 
fusion system. Using this framework, we identify similarities and 
differences between NG and other relevant models of opinion 
dynamics, such as voter model, majority game and Glauber ordering. 

Next, we present the governing rule of the opinion domain evolu- 
tion, that in the late stage of dynamics, the propagation speed of the 
opinion domain boundary is proportional to its curvature. Thus, an 
opinion domain can be considered as a mean curvature flow making 
many results of the previous works applicable here 29 " 31 . Finally we 
investigate the impact of committed agents. The critical fraction of 
committed agents found in the case of the NG on complete graph is 
also present here. We discuss the dependence of the consensus time 
on the system size, the committed fraction and the average degree. 

Results 

We begin with the definitions of the essential concepts of our model. 

Random geometric graph consists of N agents randomly distrib- 
uted in a unit square D = [0, l) 2 . Each agent has an interaction range 
defined by B r (x, y) = {(x' , y')\0 < \\(x' — x,y' — y)\\ <r},whereris 
the local interaction radius. Two agents are connected if they fall in 
each others interaction range. The choice of network topology, 
denoted as D, impacts the boundary conditions. Some studies, like 20 , 
choose the natural topology of the unit square which leads to the free 
boundary condition. In this paper, we assume that D is a torus, 
imposing the periodic boundary condition. Consequently, the opin- 
ion dynamics is free of boundary effects until the correlation length 
of the opinions grows comparable to the length scale of D. 

Microstate of a network is given by a spin vector S= (si, ■ . . ,s#) 
where s, represents the opinion of the t* individual. In the NG, the 
spin value is assigned as follows: 

{1 for A 
0 for AB . (1) 
- 1 f or B 

The evolution of microstate is given by spin updating rules: at each 



time step, two neighboring agents, a speaker i and a listener j are 
randomly selected, only the listener's state is changed (LO-NG). The 
word sent by the speaker i is represented by c, c = 1 if the word is A 
and c = — 1 if the word is B. c is a random variable depending on s ( . 
The updating rule of the NG can be written as: 

Sj(r+1) = max{min{s ; -(f)+c,l}, — l}. (2) 

Macrostate is given by n A (x, y), n B (x, y) and n AB (x, y), the concen- 
trations of agents at the location (x, y) with opinion A, B and AB, 
respectively, that satisfy the normalization condition n A (x,y) + n B (x, 
y) + n AB( x >y) = 1- We define s(x,y) = n A (x,y) — n B (x, y) as the local 
order parameter (analogous to "magnetization"), and f.i(x, y) as the 
local mean field 

l*(.X,y) = — j / / s(x\y')dx'dy'. (3) 
™ l J J B r (x,y) 

Finally, f(x,y) = - (fi(x,y) + 1) denotes the probability for an agent 

to receive a word A if it is located at (x, y). 

Through the geographic coarsening approach discussed in more 
detail in Methods, we obtain the mean-field equation describing the 
evolution of macrostate 

l^ = G ^=(l-/)-(l-/ {)* W 
while the macrostate itself is defined as 

f n A {x,y)\ 

n(X,y)=[ , , • (5) 
V n B (x,y) J 



Spatial coarsening. There are two characteristic length scales in this 
system, one is the system size (which is set to 1), the other is the local 
interaction radius r. So regarding the correlation length or the typical 
scale of opinion domains /, the dynamics can be divided into two 
stages: (1 ) Z is smaller or comparable to r; (2) l~S>r. In the second stage, 
the consensus is achieved when / grows up to 1. Figure 1 present 
snapshots of solution of the mean-field equation. They illustrate the 
formation of opinion domains and the coarsening of the spatial 
structure. 

To study the spatial coarsening, we consider the pair correlation 
function C(L, t) defined by the conditional expectation of spin cor- 
relation. 

C(L,t)=E[s(x,y,t)s(x',y',t)\ 

\\{x-x', y -y')\\=L}. 

Figure 1 implies there exists a single characteristic length scale l(t) so 
that the pair correlation function has a scaling form C(L,t) = 
C(Z = L//(f)), where the scaling function C(L) does not depend 
on time explicitly. For coarsening in most systems with non-con- 
served order parameter such as the opinion dynamics on a d-dimen- 
sional lattice, the characteristic length scale is l( t) = P". According to 
the numerical results in Fig. 2, the length scale for opinion dynamics 
on RGG at the early stage (t = 30,50) is also P, but at the late stage 
(t = 100,200,400), the length scale l(t) = p/]nt fits more precisely 
simulation results than the previous one. 

Stationary solution. Here, we find all the possible stationary solu- 
tions of the mean-field equation Eq. (4). Taking —n(x,y) = 
G(n* f) =0, we obtain 



SCIENTIFIC REPORTS | 4:5568 | DOI: 10.1038/srep05568 



2 



t=25 t=50 




Figure 1 | Snapshots of numerical solution of the mean-field equation as defined by Eq. (4). Snapshots are taken at t = 25, 50, 100, 300, the scale of 
opinion domains are much bigger than r = 0.01. Black stands for opinion A, white stands for opinion B and gray stands for the coexistence of two types of 
opinions. The consensus is achieved at t — 10". 
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7{x,y) = 



f-f+i 

V 2 -/+i 



(7) 



The eigenvalues of the linear dynamical system a = — 1 + \Jf(\ — /) 
are both negative, so n* is stable. Applying the definition of s(x, y) and 
\i(x, y), we have 

s(x,y) = n A (x,y) - n B (x,y) 

= (^(x,y))/(n 2 (x,y) + 3) 



(8) 



H(x. 



4/i(x',y; 



,>>)H 2 {x',y') + 3 



dx'dy' 



(9) 



Once we solve the above integral equation, we can retrieve the 
stationary macrostate n*(x, y) by Eq. (7). Taking fi(x, y) as a 
constant, we find three solutions fi(x, y) = ± 1 or 0. fi(x, y) = ±1 
are both asymptotically stable, while pi(x,y) = 0 is unstable. Another 
class of solution is obtained by taking fi(x, y) = j.i(x) (or similarly /((x, 
y) = n(y)). The solution consists of an even number of stripe-like 
opinion domains demarcated by two types of straight intermediate 
layers parallel to one side of the unit square D as shown in Fig. 3(b). 
With the boundary condition oo) = — 1 ,/;( + ;») = lor vice versa, 

we solve the two types of intermediate layers \x ( ^ 



as shown in 

Fig. 3 (a). The intermediate layers are of the scale r and can be placed 
at arbitrary x 0 e [0, 1). This class of solution is neutrally stable. 
Finally, there is another class of solution shown in Fig. 3(c) with 
intermediate layers both in x and y directions and opinion 
domains assigned as a checker board. This type of solution is 
unstable at the intersections of two types of intermediate layers. 
The latter two classes of solutions can be easily generalized to the 
cases when the intermediate layers are not parallel to x or y axis. Later 
we will show that in stationary solutions all the curvature of the 
opinion domain boundary has to be 0, so the solutions mentioned 
above are the only possible stationary solutions. 

In conclusion, considering the stability, the final state of the 
macrostate dynamics can be: (1) all A or all B consensus states which 
are both asymptotically stable, (2) stripe-like solution. The probabil- 
ity for the dynamics stuck in the stripe-like state before achieving full 



consensus is roughly - in analogy to similar cases in continuum 
percolation and spin dynamics 2930 . 

Two-time-scale separation of mean-field equation and compa- 
rison with other models. One important observation regarding 
the macrostate dynamics is that the change of local mean field is 
usually much slower than the convergence of local macrostate n to 
its local equilibrium n* . Let yia = n* A + 8ha, nB = n* B + 8riB, 
8n = (8nA,8riB) T . The following equation shows the exponential 
rate of convergence of n: 

— 8n(x,y) = — ( J r , \ 

St v ,y> \i-f(x,y) 

The largest eigenvalue is a — — 1 + — /) < 



{^ y) )8n(x,y). 



(10) 



So the typical 



time scale z„ of the convergence of the local macrostate is — — < 2 

which is independent of time and system size. The typical time scale 
iy of the change of the local mean field is inversely proportional to the 
propagation speed v of opinion domain boundaries, and as we will 

show later, is of the order O wri ere _R is the curvature of the 

opinion domain boundary. Therefore, zj^>z n for both long time (_R 
grows to infinity along with the time t) and big systems (in the sense 
that r<l). 

Fig. 4 shows the significant two-time-scale separation observed in 
numerical results. The equilibrium value of the local order para- 
meter, s*, can be predicted by the local mean field, 
s* — n* A — n* B = 4j.i j (/.1 2 + 3) . In Fig. 4, we present the empirical local 
order parameter s for different local mean field values /( and show 
that it is very close to its local equilibrium s*. 

odd 
Since — s(x,y) = —riA\x,y) — —nB(x,y), we have 



8t 



at 



ct 



1 



(11) 



This ODE is quite relevant to reaction-diffusion systems. On the 
right hand side, the coefficient - is easy to get rid of by scaling the 
time x = t/2. After the scaling, the first term is diffusive since (/( — s) 






(x-X Q )/r 

Figure 3 | Stationary solution, (a) Two types of intermediate layers for stationary solution /( 
slope of the intermediate layer at x = Xq is about y*lr = 1.034/r. (b) Stripe-like stationary solution, neutrally stable, (c) Checker-board-like stationary 
solution, unstable. 



-Xo 



Xq is the location of the intermediate layer. The 
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Figure 4 | Comparing the empirical local order parameter s from numerical simulation and the equilibrium s* for different local mean fields fl s. The 

solid line stands for s*. The error bars present the means and standard deviations of s that is the empirical local order parameter for the given \i in 
numerical simulations. 



is the continuum approximation of the Laplace operator on RGG 
network acting on s. The second term n AB ^i is the local reaction term. 
Though classified rigorously, it is non-local as defined in reaction- 
diffusion system, it represents a reaction in local neighborhood B r (x, 
y). The adiabaticity of the dynamics implies that the diffusion is 
much slower than the local reaction. We can obtain an approximated 
ODE for slow time scale dynamics in a closed form by estimating n AB 

by its local equilibrium 1 — n* A — n*% = ■ ^ - 



~s[x,y) = (fl-s) 



1-/(1"/) 
H(l-li 2 ) 



3 + A' 2 



(12) 



The qualitative behavior of the reaction-diffusion system is deter- 
mined by the linear stability of the reaction term 24 . In this sense, Eq. 
(12) provides clear differentiation between dynamics in our model 
and in the voter model, the majority game and the Glauber ordering. 
Taking a similar approach, we find that the voter model on RGG is 
purely diffusive, i.e. the reaction term is 0. For the Glauber ordering, 
the reaction term is tanh(/?/^) — in which /? is the inverse of 
temperature and / is the interaction intensity. Fig. 5 shows the reac- 
tion term Re((() for the voter, NG, and Glauber ordering (GO) at 
different temperatures. The majority game, NG, and Glauber order- 
ing at zero temperature have reaction terms with the same equilibria 
and stability (±1 stable, 0 unstable). Thus, the mean-field solutions 
of these models behave similarly. However, at the level of the discrete 
model, the NG on RGG will always go to a microstate corresponding 
to some stationary mean-field solution, while Glauber ordering at 
zero temperature on RGG may get stuck in one of many local minima 
of its Hamiltonian. 

Boundary evolution. The evolution of the opinion domains is 
governed by a very simple rule. The boundary of opinion domains 

propagate at the speed v that is proportional to its curvature 1/R, i.e. 

a. 

v— — . Here, a is a constant defined by the average degree (k). In 
R 



Methods, we provide a heuristic argument and using the 
perturbation method prove this relation for the mean-field 
equation, i.e. the case when (k) —> o°. This behavior is common for 
many reaction-diffusion systems and it is qualitatively the same as 
the behavior of Glauber ordering at zero temperature 23 ' 31 . 

Following the rule of boundary evolution minimizes the length of 
the domain boundary. A direct consequence of this fact is that if any 
stationary solution exists, its boundaries must be all straight (geo- 
desic), confirming our conclusion about the stationary solutions 
found in the previous paragraph. Since global topology is irrelevant 
to our derivation, this relation applies also to other two-dimensional 
manifolds. The manifold considered here is the torus embedded in 
2D Euclidean space. However, for the standard torus embedded in 
3D Euclidean space, the topology is the same but metrics are not, 
hence the geodesies are different. Therefore, there are quite different 
and more complicated stationary solutions there. Another example is 
the sphere in 3D space. On the sphere, the only inhomogeneous 
stationary solution consists of two hemispherical opinion domains, 
because the great circle is the only closed geodesic on a sphere. 

The numerical result presented in Fig. 6 confirms this relation. In 
Fig. 6, we gather 10 6 data points from numerical solutions of the 
macrostate equation, using different initial conditions, taking snap- 
shots at different times, tracking different points on the boundary 
and calculating the local curvature radius R and boundary propaga- 
tion speed v. These data points in the double-log plot are aligned well 
with the straight line with slope — 1 . The curve formed by data points 
is jiggling with some period because we implemented the numerical 
method on a square lattice, so the numerical propagation speed is 
slightly anisotropic. 

Another way to confirm this rule is to consider a round opinion 

a. 

domain with initial radius _R 0 . Given v — — , the size of this opinion 

R 2 (0) 

domain decreases as S(t) = nR 2 (t) = nR 2 , — InaX when t< — — . 
This relation is shown in Fig. 7. In simulations, the size of opinion 
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Figure 5 | Reaction term Re(//) for voter, NG and Glauber ordering (GO) at different temperatures. 
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domain S(f) is evaluated by — ^Na + -Nab J where N a ,N A b are the 

total numbers of A, AB nodes, respectively. We also observe from this 
plot that a increases with average degree (k) and converges to its 
upper bound when (k) — » °°. 

Impact of committed agents. We now consider influencing the 
consensus by committed agents. In sociological interpretation, a 
committed agent is one who keeps its opinion unchanged forever 
regardless of its interactions with other agents. The effect of 
committed agents in the NG has been studied in 15 - 13,22 . A critical 



fraction of committed agent q c = 7 — 4\/3~0.0718 is found for 
NG-LO on complete graph which is also relevant here. In our 
setting, a fraction q of agents are committed to opinion A, and all 
other agents are uncommitted and hold opinion B. The macrostate 
with committed agents is still defined by Eq. (4), but the definition of 
local mean field fi is replaced by 



fl(x,y) = ■ 



[l 9 V (*',/) + 3 



+ qdx'dy'. (13) 



Generally, q can vary on the x-y plane, but we only consider the case 
that q is a constant, i.e. the committed agents are uniformly 



10" 



> 10" 



10" 




10 



10' 

R/r 



10 



Figure 6 | Propagation speed vof the boundary of opinion domains vs. curvature radius R. Data points are gathered from 100 runs of macrostate 
equation with different initial conditions. In each run, propagation speed and curvature radius are calculated for 10000 points on the boundary. 
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Figure 7 | The size of an opinion domain S as a function of time with r = 0.05. The initial radius of the opinion domain is R 0 = 0.325. The straight line is 
for the numerical solution of mean-field equation. The dotted line and dash line are for simulation of discrete model with N = 10000 and (k) = 78.5, as 
well as with N = 20000 and (k) = 157), respectively. 



distributed. Then we reanalyze the stationary solution. Firstly there is 
a critical committed fraction q c which is exactly the same as that on 
complete graph. When q > q c , the only stationary solution is f.i(x, y) 
= 1 and it is stable. When q < q c , there are three solutions, of which 

two ji(x, y) = 1, ^\ sta ^ e ' anc ^ me tnu ~d 



fi(x,y) ■■ 



q-\ + yJq 2 -\4q+\ 



is unstable. Besides there is a class 



of stationary solutions when the committed fraction is below the 
critical. They are analogues of the stripe-like solutions in the non- 
committed agent case. The evolution of the boundary can be 
interpreted as a mean curvature flow. In such view, the fraction of 
agents committed to A opinion exert a constant pressure on the 
boundary surface from the side of A opinion domain. Similarly, 
agents holding opinion B exert a pressure from the side of B 
opinion domain. So the stationary solution will contain opinion 
domains in the form of disks with critical radius R c for which the 
pressure arising from the curvature offsets the pressure from the 

committed fraction; thus we have — ~ q. This type of stationary 

solutions are unstable, the round disk of the opinion domain will 
grow when R > R c and will shrink when R < R c . In the first case, 
when the typical length scale of opinion domains grows beyond 2R C , 
the system will achieve consensus very quickly. 

On the basis of the above stability analysis, we then analyze the 
dependence of the consensus time on the system size N, the com- 
mitted fraction q, and the average degree (k), and show our conclu- 
sions are consistent with the numerical results in Fig. 8 which for a 
given fraction a (a = 0.9 in the figure) depicts the time for a-con- 
sensus in which at least fraction a of agents hold the same opinion. 
The time to achieve a-consensus is independent of Nboth according 
to the mean field prediction and numerical plots. When q> q c the 
dynamics will converge to its unique local equilibrium fi = 1 at all 
locations simultaneously. The consensus in this case is close to that 
on the complete network, especially when (k) tends to infinity. In the 
opposite case, when q < q c , the process to consensus is twofold - 



before and after the A opinion domain achieves the critical size 2R C . 
After this criticality, the process is just the extension of the opinion 
domain driven by the mean field Eq. (4). This stage is relatively fast 
and the consensus time is dominated by the duration of the other 
stage, the one before the criticality, in which the dynamic behavior is 
a joint effect of the mean field and the random fluctuation we neglect 
in mean field analysis. Assuming the dynamics was purely driven by 
the random fluctuation, the typical length scale of opinion domains 

would have the scaling O^f 5 ^ at the early stage, hence the time scale 

of this stage would be 0(J?^) , i.e. 0(q~ 2 ). If the dynamics was purely 
driven by the mean field, the A opinion domain would never achieve 
the critical size. The actual dynamics behavior is in between the two 
extreme cases. When (k) decreases, the fluctuation level relative to 
the mean field becomes higher, leading to faster consensus. In Fig. 8, 
linear regression for the data points 0.6 < q < q c gives t c ~ q~ y in 
which y = 2.59, 2.34, 2.19 for <k > = 50, 25, 15 respectively, where 
y = 2 is the value corresponding to the purely random extreme case. 

Two observations here may have meaningful sociological 
interpretation: 

(1) When (k) — > =0, for both q < q c and q > q c , the dynamics 
behavior converges to that on the complete networks, though 
the RGG itself may be far from the complete network (with r 
kept constant, the diameter of the RGG network is much larger 
than 1). 

(2) When q < q a committed agents are more powerful in changing 
the prevailing social opinion on RGGs with low average degree. 
It is similar to the result of the previous study 22 on the social 
dynamics on sparse random networks but the "more powerful" 
is in a different sense meaning not the smaller tipping fraction 
(with longer consensus time), but the shorter time to consensus. 



Discussion 

On RGGs, the average degree of an agent (k) = nr'N is an important 
structural parameter which also strongly impacts the local dynamic 
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N=2000,<k>=15 
N=4000,<k>=15 
N=2000,<k>=25 
N=4000,<k>=25 
N=2000,<k>=50 
N=4000,<k>=50 
MF of CN 




Figure 8 | Time f c for 0.9-consensus for NG on RGG with different fractions (q) of committed agents with direct simulation on networks with average 
degrees (lc) = 15,25,50, and network sizes N = 2000, 4000. The solid black curve shows the mean field prediction of NG on the complete network(CN) as 
the limit case when (k) — » °°. The slope of the red data points near q c is —2.19. 



behavior. There are two critical values of (k): One is for the emer- 
gence of the giant component, k cl = 4.512 32 ; The other one, k c2 , only 
applicable for finite-size networks, is for the emergence of the giant 
component with all vertices belonging to it. In this paper, we only 
considered the case when (k) is above the critical value k c2 ~ In N so 
that the network is connected 26 . 

We mainly focused on analyzing the mean-field equation of the 
NG dynamics on RGG. We predicted a number of behaviors, includ- 
ing the existence of metastable states, the two-time-scale separation, 
and the dependence of the boundary propagation speed on the 
boundary curvature. We demonstrated in detail that the evolution 
of the spatial domains for the two-word LO-NG is governed by 
coarsening dynamics, similar to the broader family of generalized 
voter-like models with intermediate states 17,20,24,25 . However, there are 
still some behaviors that cannot be explained by the mean-field 
equation, such as: (i) in the large t limit, the scaling of correlation 

length is not t 1 as on the 2-d regular lattice, but P j In f; (it) the 

propagation speed increases along with the average degree (k) and 
its upper-bound is the mean field prediction, i.e., the (k) — » °° case. So 
the major limitation of the mean-field equation derived from the 
geographic coarsening approach is that it neglects the fluctuation 
among replicas (see Methods) and loses the information about (k). 
The dependence of the dynamics on (k) is left for further study. 

Methods 

Geographic coarsening approach. First, we provide the equation for the evolution of 
microstate S— (s,). Denote the probabilities for s, taking values 1, 0, —1 asp jA , p iAB 
and p jB , respectively. The master equation for spin s, is given by 



d n - 
7tP>* - 



=fi(l-piA-piB)-(l-fi)PiA 

jjPiB = (1 -/■)(! -p iA ~pm) -fipiB 



(14) 



where fi — — — is the probability for the r* agent receiving a signal A, while \i { is the 
local mean field defined as the average of the neighboring spins, 
1 

h 



E « 

'{j|(*,,»)eSr(».»)} 



(15) 



where (x h y t ) is the coordinate of the I th agent and kj is the degree of the f th agent. 
Master equations for all spins together describe the evolution of microstate. 

The motivation for geographic coarsening comes from the fact that RGG is 
embedded in a geographic space, so we may skip the level of agents and relate the 
opinion states directly to the geographic coordinates. Instead of taking into account 
the opinion of every agent, in geographic coarsening, we consider the concentration 
of agents with different opinions at a specific location. n(x,y), s(x, y),f(x, y) are 
continuously differentiable w.r.t. x,y and t We derive the equation of n(x,y) from Eq. 
(14) by taking limits, 

f <' 



-^n A {x,y) =f(x,y)(l-n A (x,y)-n B (x,y)) 
-(l-f(x,y))n A (x,y) 



— n B {x,y) = ( 1 -f(x,y) )(l-n A (x,y) - n B (x,y) ) 



(16) 



-f(x,y)n B (x,y) 



The limit is done as follows 



n A (x,y) — lim 



1 



o nc 2 NK 



E h 



(17) 



The coarsening based on the above limit is valid under either of the following two 
assumptions. The first is when the RGG is very dense (N — » °c) and c 2 N keeps 
constant. The second is when we consider K replicas of RGG with the same set of 
parameters, and the summation above is over all replicas. In addition, c 2 K keeps 
constant. Our derivation is actually based on the second assumption. However under 
the first assumption, the fluctuation is vanishing, and the dynamics behavior of a 
single run converges to the mean field result. 

Propagation speed of the domain boundary. We show here the qualitative property 

of the boundary evolution by a heuristic argument. We consider a solution with the 

form s(x,y) — g{R) where .R— \\(x — x 0 ,y— y 0 )\\ and with boundary conditions g(0) — 

— 1 andg(°o) = l,g[R) has an intermediate layer at R as shown in Fig. 3(a), so near R, 

y*(R-R) Y 
g\R)~ — . When R^$>r, using moving coordinate c, — fcx — vt where k is unit 

wave vector, x is spatial coordinate and v is the wave speed, Eq. (11) becomes 



(\ + n AB )fi-s + 2v — 5-0 . 



(18) 



Here, can be approximated by 
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f^—^J J — ( V 7 ^ 2 + r' 2 + 2Rr' cos B-R) r'dr'dd 



(19) 



yr 
6fi 



Then we make perturbation on Eq. ( 1 8) 5 — 5* + c~s, ji — + c/i, v — v* + cv and so on, 
requiring 5*(C — 0) — = 0) = v* = 5(c — 0) = 0, and obtain the equation for 

o({) 



( 1 + "Ib)/< + 'W<* + 



(20) 



At £ — 0, — 1/3, [i = fi/c, v — v/c and — s* is y*/r, the above equation becomes 
2 

-u+vy /r — 0, so 
3 ' 



2/tr _ r 



(21) 
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